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NOTATION 


r      S  Position  vectors  on  the  surface  of  the  structure 

d  f^   <Xs  Element  of  surface  area 

t  Time 

t  Time  delay 

U,  (  f^  t)  Lateral  response  of  the  structure 


Lateral  load  on  the  structure 


CD      (  r"  S^  T )      Cross-correlation  of  two  time  random  functions, 
7((r^t)and  y(S»,k)  at  points   r    and    S    respectively  (8) 

(^     (r^S^QJ       Cross  power  spectral  density  for  these  same 
functions  (9) 

Q  Angular  frequency  in  radians  per  second 

r  Frequency  in  cycles  per  second 

1  Time  length  of  a  record 

r  ^vy        Real  and  imaginary  parts  of  the  cross  power 

spectral  density  (13) 

h(f^,S,t)  Response  to  the  unit  impulse 

H(r,^^U))         Transfer  function  (Sec.    2,    Chap.   I) 

Denotes  complex  conjugate 

Hj  Discrete  representation  of  the  transfer  function 

I^  Modulus  of  the  discrete  transfer  function  (33) 
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Q  Phase  angle  of  the  discrete  transfer  function  (33) 

j  Imaginary  unit         J  =   y-l 

Oc^k),  y*(0      Discrete  notation  for    'X(r*,t)  and   y{t^t) 

Aa  Bj             Real  and  imaginary  parts  of  the  transfer  function 

d  Acceleration  of  gravity 

P  RMS  value  of  the  acceleration 

P  RMS  value  of  the  random  input  forces 

(r'(t)  Pressure  signal  used  as  pilot  signal 

N  Subscript  denoting  normalized  quantity 
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CHAPTER  I 
MATHEMATICAL  BACKGROUND 

1.     Introduction 

Although,    in  many  instances,    randonn  problems  have  reached 
a  standard  level  in  communication  engineering  and  have  been  investi- 
gated for  some  time  in  the  theory  of  turbulence  and  in  several  other 
fields,    it  is  only  recently  that  the   systematic  study  of  mechanical 
structures  under  random  loadings  has  been  undertaken  [1  through  10]. 
There  has  been  a  fair  amount  of  work  done  on  theoretical  grounds, 
most  of  it  with  the  hypothetical  white  noise  input,    but  relatively  little 
attention  has  been  given  to  the  actual  testing  of  structures  under  actual 
distributed  random  loadings   [1,  3].     The  experimental  aspect  of  the 
problem  will  be  emphasized  in  the  present  investigation.     An  attempt 
will  be  made  to  verify  the  validity  and  the  feasibility  of  a  numerical 
prediction  of  the  response  of  a  structure  to  an  actual  space-time  ran- 
dom input.     Both  the  characteristics  of  the  structure  and  the  random 
input  will  be  obtained  experimentally. 

Experimental  measurements  on  random  processes  make  in- 
tensive use  of  the  basic  concepts  of  generalized  harmonic  analysis. 


Consequently,    although  the  study  reported  here  is  experimental  in  na- 
ture,   the  nnathematical  point  of  view  of  the  problem  will  be  reviewed 
in  this  chapter  and  presented  for  a  direct  application  to  the  subsequent 
experiments. 

2.     Methods  of  Analysis 

We  shall  consider  the  space-time  response,       UL(r*,t)        ,    of  a 
linear  structure,  3  ,    subjected  to  a  lateral  space-time  load,       p(s^tj    , 
per  unit  area  distributed  over  the  entire  surface  of  the   structure.      The 
nature  of  the  load  will  not  be  specified  for  the  first  part  of  the  follow- 
ing analysis,    and  it  can  be  deterministic  as  well  as  random.      The  re- 
sponse can  be  either  a  displacement,    a  velocity  or  an  acceleration. 

Two  methods  can  be  used  to  predict  the  response.     The  first 
is  the  normal  nnode  approach  which  requires  the  knowledge  of  the 
eigenfunctions  and  the  eigenvalues  to  describe  the  characteristics  of 
the  structure,    i.e.,    nnass,    rigidity  and  damping  distribution.      This 
method  was  developed  by  several  authors  [4,  7,  8]  .      The  response  in- 
volves a  series  representation. 

A  second  method  can  be  developed  as  an  extension  to  the 
multidinnensional  case  of  the  well-known  harmonic  analysis  for  the 
one-dimensional  electrical  systems.     In  this  analysis  the  character- 
istics of  a  linear  system  are  known  either  by  the  response,     n(t)     ,    of 


this  system  to  the  unit  impulse  or  by  the  Fourier  transform  of  this 
response, 


-t-oo 


H(u))  :=    /      li(b)e'^"^dt  (1) 


The  function    H{<o)   is  called  the  system  function  or  frequency  response 
function  [11,  12,  13,  14],      The  extension  to  multidimensional  problems 
of  mechanics  has  been  investigated  by  several  authors  [6,  8].     The 
function     n(t)     now  becomes  a  function  of  two  space  coordinates,    P 
and    S  ,    nannely,         n(t"^S,t)        ,   which  is  called  response  to  the  unit 
space-time  impulse  or  sinnply  response  to  the  unit  impulse  as  before. 
Similarly,    the   systenn  function    H  (to)     becomes       n(5%S,<^)  which 

will  be  called  the  "transfer  function"   in  this  study.        In  the  two  func- 
tions,     n(r'^"s^t)  and       H(r',s',  cj)         ,    the  first  vector,  F  ,    de- 
notes the  point  of  the   structure  where  the  response  is  measured,    and 
the   second  vector,    S    ,    the  point  of  the  structure  where  the  unit  im- 
pulse is  applied.     Evidently,        n  (r^S,<4)j        is  still  the  Fourier  trans- 
form of       ntt',  S,t)  according  to  equation  (1). 

The  interesting  feature  of  the   second  method  is  that  the  struc- 
tural characteristics  are  lumped  into  a  single  function  rather  than 

^Beer  [6]   calls        \-\{f^t,cS)        the  "influence  function,"  the 
term  "transfer  function"  being  reserved  for  the  multidimensional 
space-time  Fourier  transform  of     n  (  *^,  S  ^  t ) 


being  represented  by  an  infinite  series  as  in  the  first  method.     More- 
over,   the  transfer  function       H  (r^?,  cj)      is   susceptible  to  an  experi- 
mental determination  as  will  be   shown  later.      This  is  the   second 
method  that  will  be   reviewed  here. 

3.     Determination  of  the  Response 
Let  us  assume  that  the  response  to  the  unit  innpulse,     n(r,S,  t), 
is  known.      Let     11  j    denote  the  displacement  corresponding  to  the  unit 
impulse,    n       ,    the  corresponding  velocity,    and    rl^    ,    the  acceleration. 
Then  we  have 


Taking  the  Fourier  transforms  of  (2)  and  (3),    one  readily  obtains  the 
transfer  functions  for  the  velocity  and  the  acceleration: 

H^  ~  jw  H^  (4) 

H,    .-co'^H,  (5) 

This  shows  that  once  the  response  to  the  unit  impulse  or  the  transfer 
function  is  known  for  one  type  of  response,    then  it  can  be  readily  found 
for  the  two  other  types. 


pc^t) 


Figure   1,      Load  on  a  Structure 

Let  us  now  consider  the  response      o  U  ( T^  h)     of  the  structure 
to  a  load        b(  S^h)  uS      distributed  on  aji  elemental  area   QS       ,    cen- 
tered at  point      5      as  shown  in  Figure   1.     From  the  convolution  theo- 
renn  relative  to  the  time  variable,    we  have 


du(^t)    =    /       }i(r?v)   [p(r;t.o)d?]  d 


(6) 


'•>zO 


where     •^     is  a  dummy  variable  of  integration.     The  total  response 
U  {j^j^)     due  to  a  load     D  ( f\  t)   per  unit  area  distributed  over  the  en- 
tire surface  of  the  structure  can  be  obtained  from  the  superposition 


principle,    since  the  structure  is  assumed  to  be  linear.     Integrating 
over  the  whole  surface,   we  have  * 


u(rt)  -  /  /    h(r  s,v)  p(r,  t-1?)  dsd*  (?) 


V^--o 


4.     Random  Loading 
We  shall  now  consider  a  load       b  (s'^b^     which  belongs  to  a 
random  process.      Let  us  recall  that  a  random  process  is  represented 
by  the  ensemble  of  all  possible  time  records  of  a  given  phenomenon. 
We  call  b(?,  t)      _oo<t<-j-oo         one  particular  record  of  the 

random  process.      Letting     k=  \^2.,Z,...     for  the  same     S      gives  other 
records  of  the  process.      The  ensemble, 

^  ^p(s,b)l     ,  _oo<  t<-»-oo   ,    k  =  l,2,3,... 

represents  a  randomi  process.     A  characteristic  feature  of  a  random 
process  is  that  only  statistical  properties  can  be  determined. 

We  shall  now  restrict  the  random  process  to  a  particular 
class,   namely,    a  stationary  ajid  ergodic  random  process.     Many  phys- 
ical phenomena  belong  to  this  class.     Moreover,    this  class  is  the  most 
readily  measurable.     In  a  stationary  randonn  process  the  ensemble 
averages  over  all  possible  records  are  independent  of  the  fixed  time  t_  , 


at  which  these  averages  are  taken  for  each  record.      The  ergodic  prop- 
erty means  that  a  time  average  over  one   single  record  is  equivalent  to 
the  ensemble  average.      This  last  property  is  certainly  the  most  inter- 
esting,   since  a  time  average  is  easily  realizable  electronically  and 
since  one  single  record  of  the  process  is  sufficient  to  represent  the 
entire  randonn  process.     On  the  contrary,    the  experimental  determi- 
nation of  the  ensennble  averages,    required  for  a  nnore  general  process, 
would  become  rapidly  impractical.     Verifying  that  a  randonn  process  is 
stationary  and  ergodic  is  not  an  easy  matter.      To  be  rigorous,    ensem- 
ble averages  should  be  taken  over  a  large  number  of  records  as  stated 
above  and  then  compared  to  time  averages.      This  procedure  is  often 
bypassed  and  only  partial  checks  are  nnade.     It  is  simply  verified  that 
time  averages  over  some  records  give  essentially  the  same  results. 
We  will  follow  this  procedure  here. 

5.     Correlation  Function  ajid  Power  Spectral  Density 
Both  correlation  function  and  power  spectral  density  give  a 
partial  statistical  representation  of  a  random  process.     It  should  be 
pointed  out  that  they  are  in  no  way  sufficient  to  determine  completely 
a  random  process  for  which  the  probability  density  function  also  has 
to  be  known.     However,    as  insufficient  as  they  nnight  be,    the  correla- 
tion function  and  power  spectral  density  are  useful  from  an  engineer- 
ing point  of  view,    for  which  root  mean  square  (rms)  value  and 


frequency  power  distribution  are  of  interest.     Indeed,    if  the  random 
process  is  assumed  to  have  a  Gaussian  probability  density  function, 
then  only  the   second  nnoment  is  required  for  the  determination  of  this 
function.      Therefore,    the  knowledge  of  the  cross-correlation  function, 
which  is  the  second  monnent,    will  be   sufficient  for  a  complete  statis- 
tical representation  of  a  Gaussian  random  process.     Finally,    let  us 
remark  here  that  the  output  of  a  linear  system  excited  by  a  Gaussian 
random  load  is  also  Gaussian  in  nature. 

The  correlation  function  and  power  spectral  density  can  be 
defined  independently  of  each  other  and  are  related  through  the  Fourier 
transform  pair  [11,  12] .     For  a  pair  of  time  randonn  functions,     'X  (  r^  t) 
and    _y  ( s",  t)   ,    located  in  space  at    r      and     ^     respectively,    the  cross- 
correlation  function  for  an  ergodic  process  is  defined  by 


.+  T 


cp^^Cf'.T.T)      3    l|m^    _l/       -x(rt)  i'(^t-.T)dt 


(8) 


where    T     is  a  time  delay.      The  cross-power  spectral  density,    also  for 
an  ergodic  process,    is  defined  by 


(9) 


where    to     is  the  angular  frequency. 


Between   <£>    and    m    we  have  the  Fourier  transform  relations: 


-+00 


-jCOT 


§«.(^^^)=       T.>(^'^''^)^      ^"^ 


(10) 


fx/^'^''^)   -'2^  $./^^'^)^'       d'*^ 


(11) 


Since  in  the  subsequent  auialysis  the  functions    OC    and    y    will  take  dis- 
crete values,    we  shall  introduce  the  following  notation: 


cp^^(^rT)  _^  cf;*(r) 


'l^.y('%^''^) 


C  (-) 


(12) 


r,s  =  1,2,3,... 

In  (9)  through  (12)  the  subscripts  indicate  which  functions  are  consid- 
ered in  the  correlation  or  the  spectra.      The  vectors     r     and     S     indi- 
cate where  the  functions    7<    and  y  ,    respectively,    are  located  on  the 
surface  of  the   structure.      The  superscripts    ^    and    s    are  used  for  the 
case  of  a  discrete  loading  at     I"     and    S    .      The  autocorrelation  is  ob- 
tained by  setting    -X=y   and   r  =  "^    or   ("  =  &    .     We  recall  that     (S)      (r'.s'^T) 
is  real,    but  that        ^tisV>^>^)°^      "t       ("^^     ^^  complex,    i.e., 


Cm=  '^."H^iGr.H 


(13) 


written  in  the  discrete  notation. 
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6.     Response  to  a  Random  Loading 
We  shall  apply  formula  (9)  to  the  responses    U(f^^t)  and  U(r^,k) 
at  two  different  points,    T,      and     r^    ,    viz., 


r-kX 


.-i-T 


j<^t,, 


uiroe-^^^at.  /    u(f^,ge-^"'\:it. 


(14) 


Before  using  equation  (7),    some  remarks  will  be  made   [11,  13].     Be- 
cause of  the  random  nature  of  the  loading  and  the  response,    the  upper 
limit  of  integration  for  the  variable     ^     in  equation  (7)  will  be  extended 
to  +  00    ,    since  initial  conditions  and  the  corresponding  transient  are 
not  considered  in  a  random  process.      The  load  is  assumed  to  have 
been  acting  for  a  sufficient  time   such  that  the  transient  part  of  the  re- 
sponse has  died  out.     The  lower  limit  of  (7)  can  be  extended  to  -  oo    , 
since  the  response  to  the  unit  impulse      rl(f\s*,t)     vanishes  for  negative 
values  of  time    t   .     Equation  (7)  is  now 


U(r;t)    . 


+  00 


h(ft^^)  p(s,t-»)  drdv 


(15) 


S-'V=-oo 


which  is  nnore  synnmetrical  for  the  Fourier  analysis. 

Substituting  equation  (15)  into  equation  (14),    and  rearranging 
terms,    we  have 


11 


.  +  00      r-^oo 


$uu(n.f;.'^)  = 


s/s/i),=  -«>^'>a  =  -«' 


Wf:,^>MM'J,^.^z) 


(16) 


lim  _L 


♦+T 


'+T 


j'jt, 


J«t, 


p(rt,-.,)e^    Mt,      p(^,t,-.,)e-"'^Mt 


/-T 


.dr,d$^d9d9. 


We  now  make  the  substitution 


t.=    ^   +  '^.     1 


t,  z    61+  \ 


(17) 


Using  an  abridged  notation  for  D     and  h  ,    we  obtain 


$uu('^A.'^)  = 


^rK 


+  00 


h,  e'"''  dv. 


'^.  =  -< 


«  +  o© 


-i"->a 


h,e-'     ^d9. 


'^,  =  -<x, 


lim     J_ 


•+T-V. 


hT-'Jo 


f),  e         ds-, 


y^=-T-^. 


'6,r-T-^^ 


(18) 


►  drds. 
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The  curled  bracket  is  simply  the  definition  of  the  cross-power  spectral 
density  for  the  loads      b(r^^k]  and     b(f^^b)  at  the  two  points    T,      and    ?^     , 
The  second  square  bracket  is  the  Fourier  transform  of  the  response  to 
the  unit  impulse,    i.e.,    the  transfer  function       H(r^^^^co].      The  first 
square  bracket  is  the  complex  conjugate  of  the  transfer  function 

N 

H  (f^^s^   <^)         •     Rewriting  (18),   we  have 


^aui^.Z.^)-         /     HCr,?.'-)     H(F;,?^^co)     c|)pp(2:,S,,cJ)ds;ds,  (19) 


Equation  (19)  shows  that  in  order  to  determine  the  cross-power  spec- 
trum for  the  response  <p  (njl/*^)  ^^  two  points  r,  and  r^  ,  one 
needs  to  know  the  cross-power  spectrum  for  the  load  ^  (^i>^a;^} 
at  two  points  s[  and  S^  ,  and  the  transfer  function  H  (f  ^s  ,  o)  and 
its  complex  conjugate  n(r*^'s,u))  .  At  this  time  it  is  worthwhile  to 
compare  (19)  to  the  well-known  case  of  a  one-dimensional  system  for 
which  we  have 


rr 


z 


¥jj^)  =  h^'n    (I^'^co))  (20) 

In  (20)  all  quantities  are  real  and  are  not  functions  of  space  so  there  is 
no  integration.     In  (19)  the  cross-power  spectra  and  the  transfer 
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function  are  complex  and  are   simultaneously  functions  of  two  space 
coordinates    F    and     T   .      They  have  to  be  deternnined  for  all  possible 
pairs  of  points  over  the  whole  surface  of  the  structure.      This  makes 
the  problem  much  more  involved.     We  will  now  turn  to  the  case  of  a 
discrete  loading  which  can  be  more  easily  investigated  experimentally. 

7.     Discrete  Point  Loading 
Let  us  assume  that  the   structure  is  loaded  at    n    discrete 
points      ?        (v=   1,2,3^  ...^n.^    by  concentrated  forces     j^- (t).      The  load 
per  unit  area  can  be  written  as 


p(s,t)  =  \    p^^cr-ro 


(21) 


where  oTs  -  S^j  represents  the  Dirac  delta  function  in  two-dimensional 
space.     Substituting  (21)  into  (19)  where     5,    and     S^  have  been  replaced 
by    6"^    and   ^    for  convenience,    we  have 

$uuC^/-^}  =  ^      /    /  H(rT^?:,u))   H(?,C,a))   $pp(6;e^,u)} 

^^'        s/^^  (22) 

§'(i;-SjS'(r^-i;)d6;d?^ 

Then,    since  for  a  function      X  (^) 

|(?)  S"cr-sOd?  =  _f  (rj  (23) 
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w 


e  finally  obtain,    using  the  discrete  notations  for  the  cross- spectra, 


«r 


C(-)->    ^s/-)  h:;h  ^;';h 


\'^ 


(24) 


r.^r^  ^  i,a,3,...,m. 


&.^S^.^2,3....,h. 


where     H    (u)    is  the  discrete  notation  for     H(f\S^w)    .     Also,    it  has 
been  assumed  that  the  response  is  nneasured  at  tn  locations  that  may  be 
different  from  the  number  of    n    loading  points. 

Equation  (24)  caua  be  put  into  a  matrix  form  which  is  better 
adapted  to  high-speed  computing  machines.     We  write 

H    z      H    (u)j       :     mxn   connplex  matrix  of  the  transfer  function 


H 


--  [h:M  ]  : 


conjugate  transpose  of      H 


$p=  C'>) 


:     square  complex  matrix     of  order    n    of  the 
cross-power  spectra  of  the  load  at  n    points 


m    —  ^       (wH  :     square  complex  matrix-^  of  order    lYl    of  the 


cross-power  spectra  of  the  response  at    in 


points. 


This  matrix  is  Hermitian. 


We  have 


Finally, 
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($>    H      = 


^    s. 


(25) 


§.  -   H  $  H 


(26) 


One  can  also  separate  real  and  imaginary  parts  of  the  ma- 
trices and  write 


i'u 

=  fL .  1  G. 

*r 

=  f.-jS 

H 

=    A    +iB 

in 

=  A'-iB' 

(27) 


where  the  prime  denotes  the  transpose  of  a  matrix.      Then 

*u  =  [a'-jB'][fp.jgJ[a^jb] 


(28) 


which  gives 


$^  X    a'  Fp  A  +  B'  Ff.  B  -  A'  Gp  B  +  B'  Gp  A 

4.J  (  a'g^a  +  b'g^b) 


(29) 
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8.      Particular  Case  of  Two-Point   Loading 
For  the  experimental  investigation  we  shall  consider  the  case 
of  correlated  loadings  at  two  points,    1  and  2,    on  a  structure  and  the 
response  at  one  point,    0,    on  this  same  structure.     In  this  case    ^^   is 
a  2x2  matrix,    H     is  a  column  matrix  and    ^^    reduces  to  one  real  term. 
This  means  that  the  imaginary  part  of  (29)  must  vanish. 
We  have  the  following: 


A  z 


a; 
a: 


F 


12. 


r      F 


B  = 


B 
B 


G,- 


O 


.i8' 


-g;^    o 


(30) 


Then  the  power  spectrum  of  the  response  has  the  form 


$;^Cc.)  -_  f;'h  t^s^)   +  f"(u3)  t,,(u,) 


+  F>)  TJo^)    4-  G'Cc.)  i;,(^) 


(31) 


w 


here  the       1    s  have  been  referred  to  as  transfer  coefficients  and  are 


given  by 
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T.  ^ 


-T,,= 


-  T   - 


-T.  = 
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(a:)  .  (  b:)  .  h; 

(a:)%(b:;.|h; 
a(A:A:.B:B:) 

2(a:b;-a:b:) 


(32) 


The  first  two  terms  in  (31)  represent  the  uncorrelated  part  of  the  ran- 
dom loading  and  the  last  two  terms  the  effect  of  the  cross-correlation 

of  the  same  loading. 

r  r 

If  the  modulus     M      and  the  phase  angle     9       are  used  in  the 

representation  of  the  complex  trajisfer  functions,    we  have 


h;  =  a:^jb;  ^  m; (cos e;^j sine;) 


(33) 


ajid 


T,  = 

(m:) 

K^-- 

1,.  = 

em:m;c«(6:_9:) 

T,,-- 

2m;m;  s,n(9:-9:) 

(34) 
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9.     Determination  of  the  Transfer  Function 
Before  turning  to  the  experimental  aspect  of  this  work,    we 
shall  consider  the  important  conclusion  that  can  be  drawn  fronn  the 
input-output  cross-correlation  of  a  discrete  load   p    applied  to  a  linear 
system  with  the   response    Q  .     It  is  known  that  the  measurement  of  the 
cross-correlation  of  the  input  of  a  one-dimensional  electrical  system 
with  its  output  leads  to  a  formulation  that  can  be  used  to  determine  the 
transfer  function  [11,  13].      This  result  can  be  extended  to  the  case  of 
a  multidimensional  mechanical  structure   [6],     If  a  structure  is  loaded 
at  one  single  discrete  point     S    by  a  random  load, 

p(R^h)  =   ^(r.r)  pa)  (35) 

and  the  response      U  =  U  (r^t)     is  measured  at  another  point    r     ,    then 
the  structure  caji  be  considered  as  the  usual  "black  box"  of  electrical 
engineering  and  we  have 

In  (36)      Cp     (oj)     is  the  cross-power  spectral  density  between  the  re- 
sponse   U    at  point     F    and  the  load    p   at  point      &    ,    and     •|'bp(.'^)        is  the 
power  spectral  density  of  the  load  at  point     s    .      The  first  term  is 
complex,    but  the  last  one  is  real.      Written  in  terms  of  real  and  imagi- 
nary  parts,    we  have  for  the  transfer  function    Hj(oij: 
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It  should  be  pointed  out  that  when  the  transfer  function  is 
needed  for  the  con^putation  of  the  response,    as  in  the  case  of  a  con- 
tinuous  distributed  load,    the  discrete  transfer  function    Ht(i*>)    has  to 
be  deternnined  for  a  large  number  of  locations    f^  S       ,    in  order  to  have 
an  adequate  substitute  for  the  continuous  function      H  (^r^S^iOJ.     Al- 
though it  requires  numerous  measurements  and  numerical  transforma- 
tions,   the  right-hand  side  of  (37)  can  be  determined  experimentally  as 
shown  later. 


CHAPTER  n 
EXPERIMENTAL  INVESTIGATION 

1.     Scope  of  the  Study 

The   scope  of  the  experimental  investigation  will  be  as  follows. 
First  we  desire  to  check  the  validity  of  formula  (19)  for  the  prediction 
of  the  response  of  a  structure  subjected  to  a  random  loading.      We  will 
consider  only  the  sinnple  case  of  a  random  loading  applied  simultane- 
ously at  two  discrete  points  of  the  structure  and  the  response  at  a  third 
point.      Therefore,    equation  (31)  will  be  substituted  for  equation  (19), 
thereby  replacing  the  integration  by  a  summation.     An  experimental 
determination  of  the  transfer  function  will  be  made  according  to  equa- 
tion (37),    using  a  random  loading  at  a  single  point.      Then  experiments 
with  the  two-point  loading  will  be  conducted  and  the  response  meas- 
ured.    Finally,    a  connparison  will  be  made  between  the  computed  re- 
sponse and  the  measured  response. 

Another  interesting  aspect  that  will  be  investigated  here  is 
the  effect  of  the  correlation  of  sinnultaneous  loadings  on  the  response. 
Most  existing  work  to  be  found  in  the  literature  has  been  concerned 
with  input  loadings  restricted  to  the  hypothetical  uncorrelated  white 
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noise.     In  this  case  the  imaginary  part  of  the  cross-power  spectra 
vanishes,    and  equations  (19)  and  (31)  are  greatly  sinnplified.     Actual 
random  loads  having  various  correlations  leading  to  input  power  spec- 
tra that  retain  the  imaginary  part  will  be  investigated. 

It  has  been  found  in  the  course  of  this  investigation  that  the 
conducting  of  the  experiments  and  the  processing  of  data  require  the 
setting  of  numerous  electronic  and  electronnechanical  components.     It 
was  possible  to  use   standard  electronic  equipment  for  nnost  of  the 
components,    or  at  least  to  use  this  equipnnent  with  only  minor  modi- 
fications.    On  the  other  hand,    several  electrical  circuits  and  electro- 
nnechanical  devices  were  specially  built.     Matching  the  different  pieces 
of  equipment  and  keeping  thenn  running  sinnultaneously  was  sometimes 
a  delicate  problemi.     Nevertheless,    the  experiments  proved  to  be  fea- 
sible.    Figure  2  shows  a  picture  of  the  experinnental  arrangement. 

2.     Experimental  Apparatus 

2.  1     The  Model 

The  model  tested  was  a  cantilever  beam,    7  in.    long,    1  in. 
wide  and  1/2  in.    thick,    specially  constructed  in  order  to  have  a  large 
damping  characteristic.      The  model  was  nnade  of  a  1 /4  in.    thick  plexi- 
glass core,    on  both  sides  of  which  a  1/8  in,    layer  of  a  plastic  base 
damping  material  was  cemented,    as  shown  in  Figure   3.      The  beam 
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weighs  90  grams  and  was  clamped  on  a  90-lb  stand  that  provides  a 
virtually  rigid  immovable  reference  as  shown  in  Figure  4.     Figure  5 
shows  the  location  of  the  applied  forces  at  points   1  and  2,    situated  re- 
spectively at   1   1/2  in.    and  4   1/2  in.   from  the  free  end  of  the  beam. 
The  response  was  monitored  at  point  0,    situated  1/2  in.   from  the  free 
end. 

2.  2     Loading  and  Measurement  System 

Figure  6  is  a  block  diagram  of  the  loading  and  the  nneasure- 
ments  on  the  beam.      The  two-channel  tape  recorder  no.    1  provides  an 
input  to  the  beam  of  two  random  signals  of  various  correlations.     These 
signals  were  obtained  in  the  course  of  previous  work  from  pressure 
measurements  at  the  boundary  of  an  air  jet   [15].     The  tape  was  played 
back  at  a  speed  of  15  ips,    which  is  half  the  recording  speed,    in  order 
to  n^eet  the  frequency  range  of  power  annplifier  no.    2. 

The  signal  from  channel  no.    1  was  fed  into  one  power  ampli- 
fier,   and  that  from  channel  no.    2  into  a  separate  power  amplifier. 
Amplifier  no.    1  was  able  to  deliver   1  amp  rms  into  a  15-ohm  resis- 
tance for  a  frequency  band  rajiging  from   10  cps  to  5000  cps.     Amplifier 
no.    2  was  able  to  deliver   1.5  amps  rms  into  a  50-ohm  resistance  for  a 
frequency  band  ranging  from   10  cps  to   1500  cps.      These  amplifiers 
were   standard  power  amplifiers  for  sinusoidal  waves  and  had  no  pro- 
vision for  random  wave  equalization. 
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The  output  circuit  of  each  amplifier  was  composed  of  a  load 
resistance    R   ,    a  coil    C  ,    the  resistance  of  which  was  4  ohms,    and 
a  measuring  resistance    t   of  2.  5  ohms,    all  coupled  in  series.      The 
load  resistance    R    was  added  in  order  to  reduce  the  feedback  effect 
from  the  beam  to  the  amplifier.      This  procedure  realized  a  partial 
equalization  [l]. 

Force  measurement.  --It  is   shown  in  Section  2.  3,    equation 
(38),    that  the  force  acting  on  a  coil  is  proportional  to  the  current,   1    , 
flowing  along  the  wire  of  the  coil.      Let    <-,    represent  the  current  in  the 
output  circuit  of  amplifier  no.    1.     All  the  components  of  this  circuit 
are  in  series;    therefore,    the   same  current     t    flows  through  the  coil 
and  the  resistance    *L     .      The  voltage    V[    across  the  resistance     *L,    is 
proportional  to  the  current     t,    (Ohm's  law,    Uz^i).      Consequently,    the 
voltage    V|    can  be  calibrated  to  represent  the  force       b  (tj  applied  at 
point   1.     Also,    the  voltage    V^    across    \   can  represent  the  force    P,  C^) 
at  point  2.      The  voltages   V^     and   IT     were  recorded  simultaneously  on 
tape  recorder  no.    2,    using  FM  recording  at  a  speed  of  30  ips.      This 
tape  recorder  had  seven  channels,    but  only  two  channels  were  equipped 
with  an  FM  recording  device. 

Measurement  of  the  re  sponse. --The  monitored  response  was 
the  acceleration.  The  sensor  was  a  6-gram  crystal-type  accelerome- 
ter  with  a  frequency  range  of  20  to  80,  000  cps  and  a  sensitivity  of 
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0.  114  g/mV  when  connected  across  the  preamplifier.     After  amplifi- 
cation the  acceleration  signal  was  recorded  on  tape  recorder  no.    2 
simultaneously  with  the  force  signal  no.    1.      Two  successive  runs 
were  necessary  to  record  simultaneously  the  two  force  signals  of  the 
previous  paragraph  and  then  the  acceleration  signal  and  the  force  sig- 
nal no.    1. 

2.  3     The  Coil-Magnet  System 

The  two  concentrated  forces   b  (k)  and  p„(t)  were  applied  to 
the  beann  by  means  of  a  specially  built  electrodynamic  device.      The 
basic  principle  of  this  device  is  the  one  used  in  any  comnnercial  shaJcer. 
According  to  Laplace's  law,    an  electric  current  i(k)  flowing  through  a 
helicoidal  wire,    the  coil,    placed  into  a  radial  magnetic  field  ft  ,    cre- 
ates a  resultant  force    F"(t)  acting  along  the  auxis  of  the  coil  (Fig.    7). 
This  force    F(b)  is  given  by 

F  =    (B^  L    =   k  I  (38) 


where 


F  =    F(t) 

Ui   =     induction  of  the  magnetic  field 
■^     =     length  of  the  wire 
t     =    t(t)    =     current 
K     =  VDv      =     constant 
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Figure  7.     Coil-Magnet  Systenn 
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The  main  difference  between  the  device  built  here  and  the 
commercial  shaker  is  in  the  attachment  of  the  coil  and  its  weight.     In 
the  commercial  shaker  the  coil  is  suspended  pernnanently  to  a  frame 
by  a  flexible  fixture.      The  weight  of  the  entire  moving  part  is  usually 
substantial  and  when  the  part  is  attached  to  the  structure,    it  introduces 
additional  mass  and  rigidity  that  would  have  been  prohibitive  in  the 
case  of  the  relatively  light  beam  used  here.     For  this  reason,    ex- 
tremely light  coils  weighing  only  4  grams  each  were  built  by  winding 
two  layers  of  a  0.009-in.    copper  wire  on  a  light  cylinder  of  paper,    1  in. 
in  diameter.      The  wire  was  cemented  to  the  paper  with  a  commercial 
epoxy,    and  a  conical  cap  was  attached  to  the  cylinder  to  allow  its  at- 
tachment to  the  structure.     Once  fixed  to  the  beam,    the  coils  were  an 
integral  part  of  it  (Fig.    3). 

A  radial  magnetic  field  produced  by  a  permanent  annular 
miagnet  and  its  casing  was  placed  around  each  coil.      The  gap  in  which 
the  coil  was  moving  was  approximately  1/16  in.     Each  magnet  was 
rigidly  fixed  on  the   stand  and  there  was  no  mechanical  connection  be- 
tween the  coil  and  the  magnet.      The  beam  itself  procured  the  proper 
centering  of  the  coil  playing  the   role  of  the  flexible  fixture  of  the 
shaker.      Therefore,    the  disturbance  was  reduced  to  the  weight  of  the 
coils. 
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Z.  4     Calibration  of  the  Coil-Magnet  System 

The  constant   k    of  formula  (38)  was  determined  directly  by 
a  static  test.      The  coil  was  attached  to  the  end  of  a  very  flexible  can- 
tilever beam  and  was  moving  into  the  magnet.      The  equilibrium  of  the 
beam  was  noted.     Dead  weights  were  placed  on  the  beam  vertically 
above  the  coil,    and  the  equilibrium  was  restored  by  adjusting  a  direct 
current  through  the  coil.     Figure  8  indicates  the  results  of  these 
experiments. 

In  order  to  check  the  validity  of  this  calibration  for  the  case 
when  the  current  through  the  coil  was  not  DC,    but  was  rapidly  varying 
with  time,    an  impedance  test  was  performed.     A  stiff  cantilever  beam 
with  a  large  damping  coefficient  was  built  and  driven  through  the  coil 
by  a  sinewave  at  the  resonance  frequency  (T    =  500  cps)  of  the  funda- 
mental flexion  mode.     It  is  known  that  near  the  resonance  one  can  con- 
sider a  structure  as  a  one  degree-of-freedom  system.      The  relation 
between  the  complex  driving  force     Jf    and  the  complex  response    A-     at 
the  point  of  application  of  the  force  is  given  by 

[(-  |A  co^  ->-  Y)  -t-  j   b  '^1'X  -  ^  (39) 

where  |jI  is  the  generalized  mass,  TiT  the  generalized  rigidity  and  D 
the  generalized  damping  constant.  At  the  phase  resonance  frequency 
W^      ,    we  have 
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Figure  8.      Calibration  of  Coil-Magnet  System  by  Static  Test 
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\XU)^    --f  (40) 


If  ^    is  the  complex  acceleration  and   W    the  reduced  damping  coeffi- 
cient,  we  have 


-LO'-X     ^    A  0(  = 


(41) 


Therefore,    equation  (39)  becomes  (shifting  fronn  complex  to  real 
representation) 

J>     =  2  jJl  A'*  (42) 

where    J*     and  Jt     are  either  the  nnaximum  value  or  the  rms  value  of 
the  force  and  the  acceleration.     Formula  (4Z),    together  with  (38),    was 
used  to  compute  the  acceleration  jz  ,    once  the  force    fi    >    the  damping 
coefficient    o(    and  the  generalized  mass   kX    were  measured.      This  conn- 
puted  value  of  the  acceleration   J\    was  then  compared  to  the  measured 
value  by  means  of  an  accelerometer. 

The  damping  coefficient    0(    was  determined  by  a  transient  at 
the  natural  frequency     C       : 


loq  >f„       -  loq  ■&- 
«=   0.7  34-  _f__511^: i_!: 


(43) 


where     iJ-     is  one  nnaximum  amplitude  of  the  transient  and    ■»?;_         the 
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maximum  amplitude  after    n    half- oscillations.     For  the  calibration 
beam  used  here  we  obtained    0<  =  0.  017. 

The  generalized  mass   \k    was  determined  by  the  method  of 
variation  of  frequency  [16]   resulting  from  small  additional  masses  on 
the  vibrating   structure.     If  a  small  additional  mass     A  (i   is  added  to 
the  structure  at  the  point  of  measurement  of  the  response,    the  struc- 
ture has  then  a  new  resonance  frequency    !>,-•-  ^\^     •     From  equation 
(40)  we  have 


Aii  +  2   Ml    =  O  (44) 


or 


u  =  _:^  A^!_  (45) 


The  frequency  shift    A-f„   was  measured  by  a  beating  process  with  the 
original  natural  frequency    C     before  the  additional  masses  were  added. 
The  experiment  gave    W  =  456/g  in  grams/nn-sec"    ,    where  g  is  the 
acceleration  of  gravity.     Finally,    the  relation  between  the  acceleration 
and  the  current  through  the  coil  for  the  calibration  beam  is 


A  -   }L  i  (46) 

1  15.4 
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where 

•^    =     the  acceleration  in  number  of  g' s 

<-      =     the  current  in  amps 

K     =     the  constant  given  in  Figure  8. 
Figure  9  gives  a  comparison  between  the  connputed  values 
according  to  (46)  and  the  measured  values  of  the  acceleration.     It  can 
be  seen  that  the  constant   K    measured  by  the  static  calibration  gives 
a  satisfactory  value  for  dynamic  loading,    at  least  for  the  frequency 
tested,    i.e.,    500  cps.     It  was  assumed  that  the  agreement  was  also 
acceptable  for  the  whole  range  of  frequencies  of  the  tested  miodel,    i.e., 
from   100  to  1500  cps. 

3.      Testing  Procedure 

3.  1     Relative  Sign  of  Forces  and  Acceleration 

Measuring  the  force  by  the  voltage   'U'   across  the  resistance    t 
(Fig.    6)  does  not  provide  the  sense  of  the  applied  force.     This  sense 
depends  upon  the   relative  sense  of  the  induction     ^    and  the  current     ^ 
in  the  coil,    and  not  upon  the  resistance   *l   .     It  is  important  to  measure 
the  force  and  the  corresponding     acceleration  with   the    proper    rela- 
tive sign.      This   sign  was  determined  by  sending  through  the  coil  a 
sinewave  at  a  frequency  near  the  fundamental  natural  frequency  of  the 
beam.     Force  and  acceleration  have  to  be  in  phase  below  the  natural 
frequency  and  out  of  phase  above. 


37 


Acceleration  in  g's 


15 


lo 


Computed  acceleration 
for   coil-magnet  no.    2 


leration 
gnet  no.    1 


± 


o.z 


Oi- 


O.fe 
Amps 


o.a 


10 


Figure  9.     Verification  of  Coil-Magnet  Calibration 
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3.  2    Experimental  Determination  of  the  Transfer  Function 

Only  one  force  at  a  time  was  applied  in  this  series  of  experi- 
ments.     Tape   recorder  no.    1  was  played  back  and  only  channel  no.    Z 
fed  into  amplifier  no.    2,    giving  a  random  force  at  point  2  of  the  beam. 
Force  and  acceleration  signals  were  recorded  simultajieously  on  tape 
recorder  no.    2  for   15   sec. 

Amplifier  no.    2  was  then  disconnected  from  coil  no.    2  and 
connected  to  coil  no.    1.      The  whole  procedure  was  repeated,    force  and 
acceleration  signals  being  recorded  for  points   1  and  0  on  tape  recorder 
no.    2.      The  setting  of  the  attenuation  for  the  acceleration  annplifier  and 
the  force  attenuator     f\\     (Fig.    6)  were  carefully  noted  for  subsequent 
determination  of  the   rms  values.      These  two  tests  provided  enough 
data  for  the  computation  of  the  transfer  function  corresponding  to 
points   1,    2  and  0. 

3.  3     Two-Point  Loading  Tests 

In  this  series  of  tests  both  channels  of  tape  recorder  no.    1 
and  both  amplifiers  were  used  in  order  to  deliver  simultajieously  two 
correlated  forces,      p  (l)  and    p^(t)  ,    at  two  points,    1  and  2,    of  the  beam. 
The  desired  data  to  be  recorded  were  the  two  forces    |3,(,t)    and    p^  (t) 
and  the  acceleration    Cl(b)  taken  at  point  0.     As  stated  above  in  Section 
2.  2,    two  runs  were  necessary  to  record  these  three  signals.     Figures 
10  and  11  show  a  picture  of  a  sample  of  the  random  forces  and  the 
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Figure   10.      Typical  Force  Signal  (Points  1  and  2) 


Figure   11,      Typical  Signals  of  Acceleration  (Point  0) 
and  Force  (Point  1) 
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acceleration  displayed  on  an  oscilloscope.     Attenuations  for  the  ac- 
celeration and  the  forces  were  carefully  noted  as  above  for  the  deter- 
mination of  the  rms  values. 

3.  4    Calibration  of  the  Electronic  Chains 

The  calibration  of  the  electrical  network  for  measurement  of 
the  two  forces  was  done  by  employing  a  500-cps  sinewave  rather  than 
tape  recorder  no.    1.      The  gain  of  the  amplifiers  was  adjusted  to  have 
3  volts  rms  across  each  resistance    *l,    and    h.^  .     These  two  voltages 
It       and    IT.       were  measured  with  a  standard  vacuum  tube  AC  volt- 

C,  C2, 

meter  and  recorded  on  tape  recorder  no.    2  with  a  known  attenuation. 
They  were  equivalent  to  the  following  forces  on  the  coils  and  could  be 
interpreted  as  forces: 

1^     =     0.  216  kg  on  coil  no.    1 

V^     =     0.  290  kg  on  coil  no.    2. 
Likewise,    the  calibration  of  the  electric  network  for  measurement  of 
acceleration  was  acconnplished  by  means  of  a  sinewave  at  500  cps.     An 
rms  voltage  of  6  1  mV  was  substituted  for  the  output  of  the  accelerome- 
ter.      The  output     IT^.     of  the  voltage  annplifier  was  recorded  on  tape  re- 
corder no.    2.     According  to  the  calibration  of  the  acceleronneter,    the 
voltage    "^Co     recorded  on  the  tape  is  equivalent  to  an  acceleration  of 
6.95  g's. 


41 


A  comparison  of  the  random  signals   ^,(^)    ,     Pay)  'Lnd    0.(k)  re- 
corded on  the  tape  with  the  calibration  signals    V^      ,  V^       and    iT^     pro- 
vides the  rms  value  for  these  random  signals  at  playback.     It  must  be 
noted  here  that  the  rms  of  the  random  signals  cannot  be  measured  by 
the  standard  AC  voltmeter.        A  true  squaring  device  is  necessary. 
This  was  done  by  using  the  combination  multiplier-DC  voltmeter  de- 
scribed in  the  following  section. 

Eight  different  correlated  loads  divided  into  two  series  were 
tested.      The  correlation  of  these  loads  is  shown  in  Figures  Z5  and  39 
in  Chapter  III.     The  case  of  uncorrelated  load  was  also  investigated 
(Fig.    35,    Chap.   lU). 

4.     Data  Reduction 

4.  1     Cross-Correlation  and  Cross-Power  Spectral  Density 

Direct  experimental  determination  of  the  cross-correlation 
ajid  the  cross-power  spectral  density  are  both  possible.     The  cross- 
correlation  has  the   simpler  experimental  determination.     For  an  er- 
godic  process,    the  latter  determination  requires  a  delay  line,    a  mul- 
tiplier and  an  averager.     If  the  cross-power  spectral  density  is  also 
needed,    as  it  is  for  this  case,    it  can  be  obtained  by  talking  the  Fourier 


A  standard  AC  voltnneter  is  a  rectifying  device,    the  calibra- 
tion of  which  depends  upon  the  wave  shape.     It  is  calibrated  only  for  a 
sinewave. 
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transform  of  the  cross-correlation,    which  can  be  done  numerically 
(Sec.    5).      The  amount  of  computation  becomes  very  large  when  the 
cross-correlation  curve  is  oscillating  and  approaching  zero  at  a  slow 
rate.     Only  high-speed  computing  machines  can  then  handle  the  prob- 
lem.    Electronic  systems  have  also  been  devised  for  this  purpose 
[13,  17]. 

The  direct  experinnental  determination  of  the  cross-power 
spectral  density  is  also  possible   [  18]  ,      The  method  uses  a  spectrunn 
analyzer  as  a  basic  instrument.     A  spectrum  analyzer  can  measure 
only  real  power  and  can  process  only  a  single   signal  at  a  time,    ex- 
cluding any  cross  measurement  between  two  signals  as  required  in 
the  cross-power  spectral  density.      This  difficulty  has  been  overcome 
by  using  adding  and  subtracting  devices  together  with  a  delay  line  [  18] 
This  method  seems  to  be   somewhat  miore  difficult  than  the  direct  de- 
termination of  the  cross-correlation.     If  the  cross-correlation  is 
needed,    the  inverse  Fourier  transfornn  can  be  taken  nunnerically,    but 
the  numerical  integration  is  in  no  way  quicker  than  the  one  in  the  pre- 
vious method.      The  direct  determination  of  the  cross-correlation  will 
be  performed  here. 

4.  2    Normalization 

All  the  correlation  curves  and  the  spectral  density  curves 
reported  in  this   study  have  been  normalized.      This  procedure  permits 
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a  better  comparison  between  the  different  curves.     The  normalization 
is  such  that  all  the  autocorrelation  functions  are  equal  to  unity  for  a 


zero  time  delay:    viz.  , 


-?:» - ' 


(47) 


This  means  practically  that  each  signal  has  been  divided  by  its  rms 
value,    since 


"♦T 


I  gcx     /         x_»oo   / 


(48) 


where  the  bar  denotes  mean  value.      Therefore,    the  relation  between 
normalized  cross-correlation  and  true  cross-correlation  is 


Cw 


?::  w 


^(t_T)  yft) 


^     V^fVlyct)]'       \/[^w]'V[:/ct)]' 


(49) 


where 


is  dimensionless. 


For  the  power  spectra  we  have 
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fs. 


Cm 


"   Vw^V  [«')]'' 


(50) 


where 


•■S/ 


has  the  dimension  of  tinne  . 


N 


4.  3     Length  of  the  Sannple 

The  con^putation  of  the  cross-correlation  was  done  according 
to  its  definition,    equation  (8),    and  is  repeated  here  in  its  simplest 
notation: 


f.y^^^    - 


=    1  im 


T-oo     2T 


(51) 


This  definition  requires  the  length  of  the  records     'X(l)  and   y(t)    to  be 
infinite.     For  practical  consideration,    only  truncated  samples  are 
available.      Therefore,    fornnula  (51)  is  replaced  by 


'+T 


^%,W    --^\       '^(t-X)^(lr)dt 


(52) 


where    "T   is  the  total  time  length  of  the   sannple.      Choosing  the  proper 
time  ~r  is  a  delicate  problenn.     It  was  assumed  here  that  if  the  time  T 
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was  much  larger  than  the  period  of  the  smallest  frequency  component 
of  the  truncated  record,    then  formula  (52)  was  a  good  substitute  for 
(51),     It  was  also  kept  in  mind  that  the  time  delay   T    should  not  exceed 
a  reasonable  percentage  of    T   which  has  been  suggested  to  be  about 
5  per  cent  [  19] . 

Different  tinnes   T    were  tried  for  the  determination  of  (5Z). 
The  cross-correlation  functions  were  first  nneasured  for  ZO-sec  sam- 
ples,   playing  back  and  rewinding  the  tape  for  each  chosen  value  of    T  , 
Then  subsamples  of  8,    4  and  Z  sec  were  cut  from  the  original  samples. 
With  such  shorter  samples  it  was  possible  to  make  loops  that  could  be 
played  back  continuously  on  the  tape  recorder.     Very  little  variation 
was  found  in  the  computation  of      (D     (z)     for  the  different  loop  lengths. 
The  time  delay  f    never  exceeded    T  =  20  nnillisec.     The  ratio   '^/y  is 
then  0.5  per  cent,   which  is  much  lower  than  the  value  suggested  [19]. 
It  was  then  concluded  that  a  4-sec  loop  could  be  adopted  as  a  good  rep- 
resentative sample  of  the  randonn  signals,    and  that  (5Z)  would  give  a 
good  estimate  of  the  cross-correlation  function. 

4.4    Correlator 

A  simple  analog  computer  was  assembled  for  the  computation 
of  (52).      The   sequence  of  operation  is 

1.     Introduce  a  time  delay  in  one  of  the  signals   %    or  y 
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Z.      Perform  a  multiplication 

3.     Take  a  time  average  of  the  product. 

Figure   IZ  shows  a  schematic  of  the  computing  system.     The 
two  signals    ^(t)   and  y(t)   are  picked  up  sinaultaneously  fronn  the  play- 
back of  the  FM  tape  recorder  no.    Z.     One  of  the  signals,    depending 
upon  the  sign  given  to   T  ,    is  fed  into  an  audio  delay  line  that  provides 
a  given  time  delay,    variable  by  steps  of  0.  1  millisec  from  0  to  Z  nnilli- 
sec.     Then  the  signals  are  amplified  and  the  multiplication  is  per- 
fornned  by  a  commercial  multiplier, 

^m(t)  =   -AO-T)  yCO  (53) 

Figure   13  shows  the   square  of    ?^(t)  and  Figure   14  shows  the  product 

The  third  operation  is  nothing  more  than  determining  the 
mean  value  of  the  product  signal     ni(t).     Electrically  this  means  the 
measure  of  the  DC  component  of  ni(t),   which  can  be  done  by  a  DC 
voltmeter.     However,    the  product  operation    *(Jb)    by   y(t)    introduces 
low-frequency  components  in     m(t)  ,    although  they  did  not  exist  in 

The  product  signal  ^mCt)   can  be  considered  as  having  the 
following  Fourier  series  expansion: 


T 


where     171^      is  the  mean  value  measured  by  the  DC  voltmeter. 
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Figure   13.      Typical  Squaring  of  Random  Signal 


Figure   14.      Typical  Product  of  Randonn  Signals 
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'X     and    y   ,    and  makes  the  meter  oscillate.      To  smooth  out  the  reading, 
a  resistance  capacitor  box  was  built,    the  effect  of  which  was  to  elimi- 
nate the  ripple  of  the  DC  component  by  shorting  out  all  AC  components. 
This  box  was  placed  between  the  output  of  the  multiplier  and  the  meter. 
The  loading  of  the  capacitors  of  the  box  required  about  20  sec. 

4, 5     Procedure 

The  use  of  a  loop  greatly  reduced  the  processing  time  by 
eliminating  the  rewinding  time  of  the  tape  recorder  and  having  the  RC 
box  almost  continuously  loaded.     According  to  formula  (49),    three 
averages 


a(b-T)  v/CO     ,    [xCt-r)]"    ^    [y(t)]  (54) 

are  necessary  to  obtain  the  normalized  correlation.     The  procedure  is 
to  fix  a  value  of    T   (negative    T  being  introduced  by  delaying  the  second 
signal   ^(t)  ),    then  perform  successively  the  three  averages  (54).     The 
actual  measurement  of  the  two  mean  square  values  and,    consequently, 
the  division  in  (49)  were  avoided  by  adjusting  the  gain  of  each  amplifier 
in  such  a  manner  that  the  meter  would  read  iinity 


[K(t-r)]'    .:    [yCL)]'=  I  (55) 
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in  performing  the  mean  square  operation.      The  reading  of  the  mean 
value  of  the  product     mCO  gives  directly  the  normalized  correlation 

X 


.      The  method  was  even  rendered  more  rapid  by  using 


two  multipliers  simultaneously,    one  for  the  product  and  the  other  for 
the  mean  square.      The  procedure  was  repeated  for  each  value  of    T  , 
and  the  correlation  curves  plotted  point  by  point  on  large-scale  paper. 
Constructing  a  correlation  curve  required  an  average  of  100  meter 
readings. 


Figure   15.     Recopying  the  Delayed  Signal  on  the  Loop 

The  delay  line  used  in  Figure  12  has  a  range  of  only  2  milli- 
sec,    but  the  correlation  curves  reported  are  extended  to  20  millisec 
in  some  cases.     This  long  time  delay  was  done  by  copying  the  2- 
millisec  delayed  signal    7(1-2)  and  the  undelayed  signal  v/(t)from  two 
channels  of  the  loop  to  two  other  channels  on  the  same  loop  (Fig,    15). 
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When  played  back,    these  two  new  signals  have  an  original  time  shift  of 
2  millisec  that  comes  in  addition  to  the  value  of    T   set  on  the  delay  line. 
This  procedure  can  be  repeated  at  will  to  produce  time  delay  up  to  the 
desired  value.     It  should  be  pointed  out  that  this  was  possible  only  with 
an  FM  recording  technique,    since  this  technique  erases  whatever  has 
been  recorded  on  the  track  before  a  new  signal  is  to  be  recorded. 

5.     Numerical  Computation  of  the 
Power  Spectral  Density 

All  the  power  spectral  densities  (and  cross-power  spectral 
densities)  reported  in  this  experimental  study  have  been  obtained  by 
taking  the  Fourier  transform  of  the  autocorrelation  (and  cross- 
correlation)  curves.      Let  us  rewrite  (10)  in  the  discrete  notation; 
omitting     %   and    y   for  simplicity,    we  have 


■Voo 


iuT: 


$•"'(0,)  ^  /       cp'n^)e"^    'dr 


(56) 


We  now  separate  real  and  imaginary  parts  in  (56)  and  use  the  symme- 
try property 

cf-*(r)  -  cp^'C-r)  (57) 

Using  the  notations  of  (13)  for  the  real  and  imaginary  part  of      <p     (lp)  , 
we  have 
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Figure   16.     Fourier  Transform  Integration 


F    (U)):. 


c^''(T)  +  cp'^r) 


C06  uj-c  dr 


(58) 


G'''(co)  =    /      [cp^'Ct)    -    cp^Cr)]    sin    a>r   d' 


(59) 


The  integration  was  performed  numerically  by  approxinnating 
the  curve      CD    C"^)   4-  di    Ct)    by    n    horizontal  segments  of  straight  lines 
of  length    AT     (Fig.    16).     The  interval  of  integration  was  also  split  into 
n    equal  subintervals: 
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4?  rT,^^ 
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•■      2 


F  Coo)  ~  cD^o)  I     Cos  uizdz  +...  + <V   (r.)/       cosojx  dr  -h...  (60) 

'r.   AJ 

Then  the  cosine  was  integrated  into  each  interval.     The  same  proce- 
dure was  followed  for  (59).     Finally  we  have 


■I'!', 


F   (i)  --L  sin  Trrar  , 

^  TTf  ^  1 


^'\o)  +^  [cp''*(r.)  +  c[)^'c^,^]  COS  aTTJ-r,  L     (61) 


(62) 


where  it  is  recalled  that  the  indices    9<    and    y  have  been  omitted  for 
sinnplicity.     These  formulas  derived  for  the  cross-power  spectral 
densities  are  evidently  valid  for  the  ordinary  power  spectral  density 
for  which  (62)  will  be  identically  zero.     The  numerical  values  of  <p    \X) 
were  read  from  the  large-scale  correlation  curves,    and  the  interval 
At"    was  chosen  as  0.  08  nnillisec.      The  numerical  computation  was 
done  on  an  IBM  650  computer. 


CHAPTER  m 


INTERPRETATION  OF  THE  RESULTS 


1.     Transfer  Functions 
The  transfer  functions     H|(|)  and    H^(^)    have  been  determined 
by  the  cross-correlation  technique  explained  in  Section  9,    Chapter  I, 
Autocorrelation  auid  cross-correlation  functions  were  measured  ac- 
cording to  the  procedure  of  Section  4.  4,    Chapter  II,    and  are  shown  in 
Figures   17,    18  and   19.      Then  the  power    spectral  and  cross-power 

spectral  densities,      $",,(().    <|J,p(f)       -f^pC^       '$^pU^       '^^^« 
computed  numerically  according  to  (61)  and  (62).     Finally,    the  trans- 

0  0 

fer  functions     H  (f)  and   H,(f)  were  obtained  by  using  formula  (37)    and 

are  represented  in  Figures  ZO  and  21.     The  modulus  and  phase  angle 

representation  has  been  chosen  for  ease  of  plotting.     Normalized 

values  of  the  modulus  are  plotted  in  Figure  20.     It  can  easily  be  shown 

by  using  (36)  and  (50)  that  the  actual  value  of  the  transfer  function    H    (f) 

is  related  to  the  normalized  value      H    (■(]      by 

-In 


Hlcn-  Kw 


M    p 


(62) 


w 


here    1      is  the  rms  value  of  the  acceleration     clCOi    and    P    is  the  rms 
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Figure  20.     Modulus  of  Transfer  Functions  H  (j)  and  H  (f) 
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Figure  21.     Phase  Angle  of  Transfer  Functions    H,(f)and  Hjj(-f) 
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value    of   the    applied  force     pCO  •       The    ratio      i/  P  is  indicated 

in   Figure  20. 


The  modulus     n 


and  the  angle     9        for  point   1  do  not  ex- 


hibit very  nnuch  variation  with  the  frequency  except  at  the  two  ends  of 


the  spectrum.     But  the  modulus 


h: 


for  point  2  has  a  sharp  maximum 


for  a  frequency  of  about  500  cps,    which  corresponds  to  the  natural 
frequency  of  the   second  mode  of  the  beam.      The  angle      Dj^     also  has  a 
rapid  variation  for  that  sanne  frequency.      The  values  obtained  for  the 
frequencies  above   1200  cps  nnust  be  considered  with  caution  since  the 
input  force  spectra     r_^      and      r         and  the  cross  spectra    ®  and 

a>  all  have  small  amplitudes  in  this  region.     Therefore,    equation 

(37)  is  close  to  the  undeterminate  form  0/0  which  evidently  gives  a 
poor  accuracy.     Nevertheless,    a  tendency  for  a  second  peak  at  a  fre- 
quency in  the  neighborhood  of  1400  cps  can  be  detected.      The  transfer 
coefficients,     I,,      ,    T^^    ,     Ij^    ,      \^^    ,    corresponding  to  the    Hs  have 
been  computed  according  to  (32),    and  are  represented  in  Figures  22 
and  23.     We  note  that  in  these  figures  actual  values  of  the     Is  are  re- 
ported rather  than  the  normalized  values. 

2,     Loading  of  Series  I 
Figure  24  represents  a  typical  autocorrelation  of  one  force 
for  a  two-point  loading.     Figure  25  represents  the  various  cross- 
correlations  of  the  input  forces  of  series  I,    applied  simultaneously  at 
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points   1  and  2.      These  forces  correspond  to  actual  pairs  of  air  pres- 
sures measured  at  different  locations  downstream  from  an  air  jet  [18]. 
The  effect  of  these  different  loadings  on  the  acceleration  of  the  model 
is   shown  by  the  autocorrelation  curves  in  Figure  26.     A  reasonably 
rapid  decay  of  these  curves  has  been  obtained  because  of  the  large 
damping  coefficient  originally  introduced  into  the  model.      The  Fourier 
transform  of  all  these  curves  leads  to  the  corresponding  cross  spectra. 
Each  of  the  figures,    27  through  36,    represents  the  normalized 
spectra  and  cross  spectra  for  the  input  loadings  and  the  corresponding 
accelerations.     Actual  values  of  these  spectra  and  cross  spectra  are 
obtained  according  to  formula  (50).      The   spectra  and  the  cross  spectra 
of  the  input  forces  exhibit  occasional  fluctuations  in  terms  of  the  fre- 
quency which  did  not  exist  in  the  input  pressures  from  the  tape  before 
amplification  (Fig.    5).     These  fluctuations  are  believed  to  be  due  to  the 
feedback  effect  of  the  beam,    through  the  coil,    on  the  output  of  the  power 
amplifiers,    nos.    1  and  2.     It  must  be  noted  here  that  since  the  coils    C, 
and    C„     are  each  in  series  with  the  nneasuring  resistances     't^     and 
•^-     ,    the  voltages    V,     and    V     across     1,     and    ^^    give  the  real 
forces  acting  on  the  beam,    i.e.,    the  intended  forces  minus  the  feed- 
back forces.     Therefore,    the  forces  applied  at  points  1  and  2  are  not 
exactly  the  same  as  the  pressures      JjCO  and    S^i-)    from  the  tape,   but 
the  departure  of    b(t3  and  pj^(t)from    ij(l)    and    S^(^)  is   small  enough  to 
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Figure  27.     Power  Spectral  Densities  of  Forces,    Case   1 
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Figure  28.     Power  Spectral  Densities  of  Acceleration,    Case   1 


69 


Power  spectral  density 


sec 


0 

-^1 L 


-2-22 


Curves 

rms 

2-11 

0.214  lb 

2-22 

0.311  lb 

2-12 

0.258  lb 

2-11 


500 


looo 


1500 


Frequency,    cps 
Figure  29.     Power  Spectral  Densities  of  Forces,    Case  2 
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Figure  30.     Power  Spectral  Densities  of  Acceleration,    Case  2 
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Figure  31.     Power  Spectral  Densities  of  Forces,    Case  3 
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Figure  32.     Power  Spectral  Densities  of  Acceleration,    Case  3 
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Figure   33.     Power  Spectral  Densities  of  Forces,    Case  4 


74 


0 


Power  spectral  density 


]  f 


100  t- 


10 


0.1 


10' 


sec 


Curves 

rms 

Measured 

4,  10  g 

•      Computed 

4.30  g 

500 


1000 


1500 


Frequency,    cps 
Figure  34.     Power  Spectral  Densities  of  Acceleration,    Case  4 
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Figure  35.     Power  Spectral  Densities  of  Uncorrelated  Forces,    Case  5 
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Figure  36.     Power  Spectral  Densities  of  Acceleration,    Case  5 
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produce  only  a  small  change  in  the  intended  cross -correlation.     Never- 
theless,   since  we  nneasure  the  actual  forces  applied  at  points   1  and  Z, 
we  can  forget  about  the  feedback  effect.      The  correlations  of  Figure  Z4 
and  also  of  Figure  39  represent  the  real  loading,    which  is  very  close 
to  the  intended  one. 

Figure  27  represents  the  power  spectra  and  cross-power 
spectrum  of  the  loading  denoted  as  case   1.     In  this  case  the  sanne  pres- 
sure signal    J"(t)  was  fed  simultaneously  into  the  two  power  amplifiers, 
nos.    1  and  2  (Fig.    6).      This  procedure  gives  the  maxinnunr^  cross- 
correlation  between  the  two  forces   jp^Ck)    and    pj(t)  ,    i.e.,    a  zero  imagi- 
nary part  in  the  cross-power  spectrum.     However,    due  to  the  phase 
shift  between  the  two  amplifiers  (they  were  of  different  construction) 
and  the  feedback  effect,    it  can  be  seen  from  Figure  27  that  the  imagi- 
nary part   1-12  is  not  zero.     Nevertheless,    the  correlation  is  strong, 
as  indicated  by  the  real  part  1-12.     The  spectrum  of  the  measured  ac- 
celeration at  point  0  is  represented  in  Figure  28  by  the   solid  line.      The 
plot  in  this  figure,    as  well  as  for  the  subsequent  spectra  of  the  accel- 
eration for  other  loadings,   has  been  made  on  semilog  paper  because  of 
the  sharp  variations  of  the  spectrum  in  the  neighborhood  of  5  00  cps. 
The  dots  on  Figure  28  represent  the  computed  spectrum  according  to 
formula  (31),    using  the  transfer  functions  determined  experinnentally 
above.      The  connputation  of  (31)  requires  the  knowledge  of  the  rms 
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values  of  all  the  quantities  involved  besides  the  normalized  spectra 
and  transfer  functions.     Indeed,    the  response  depends  not  only  upon 
the  nature  of  the  two  forces     b(_i)  and  p^(t)  ,    but  also  upon  their  rela- 
tive strength,    ajid  upon  the  size  and  stiffness  of  the   structure.      The 
rms  values  of  the  forces    F^  =  l^p^Ct)]^  .    P^   -^  [jD^(b)]^  '  and   P,^-\[^ 
are  indicated  in  Figure  27.     Formula  (50)  was  used  to  obtain  the  actual 
spectra  of  the  input  forces. 

O  0 

Once  the  spectrum  of  the  acceleration    ^^^C-f)  had  been  com- 
puted by  (31),    then  the  computed  rms  value  could  be  determined  by 
integrating  the  spectrum  over  the  whole  range  of  frequency,    viz., 


=\/[awr=u/  <!>::«)  df 


The  computed  rms  value    T    is  indicated  in  Figure  28,    and  the  dots 
represent  normalized  values  of  the  computed  spectrum. 

The  agreement  between  the  experimental  curve  and  the  com- 
puted values  is  not,    in  general,    excellent,    but  is  nevertheless  accept- 
able if  one  considers  the  nunnber  of  successive  experimental  nneasure- 
ments  and  numerical  computations  involved  between  the  test  runs  and 
the  curves  in  Figure  28.     Running  several  identical  tests  rather  than 
only  one  of  each  case  as  reported  here,    would  certainly  show  the  mag- 
nitude of  the  dispersion  and  improve  the  agreement  by  averaging  over 
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the  number  of  runs.      The  computed  rms  value  and  the  rms  value  de- 
termined experinnentally  differ  by  only  3  per  cent  which  is  good  but  is 
considered  as  an  exception,    since  this  difference  ranged  between  5  and 
10  per  cent  for  the  other  cases  of  loading. 

The  general  shape  of  the  power  spectrum  is  reasonably  well 
predicted,    at  least  for  the  values  of  the  frequency  below   1000  cps. 
Above  1000  cps  the  agreement  is  poor  as  can  be  expected,    since  the 
transfer  function  is  itself  not  determined  with  great  accuracy  in  this 
region. 

Figures  29  through  36  represent  different  cases  of  correlated 
loadings  that  were  tested  in  series  I  and  the  corresponding  measured 
and  computed  spectra  of  the  acceleration.     It  can  be  seen  that  all  out- 
put spectra  have  a  nnaximum  at  about  500  cps,    but  the  shape  of  the 
curves  varies  fronn  one  case  to  the  other,    showing  the  effect  of  the 
various  loadings.     Analogous  to  the  situation  in  case   1,    the  general 
shape  of  each  curve  is  reasonably  well  predicted  except  for  the  fre- 
quencies above  1200  cps.     The  computed  points  fall  off  the  measured 
spectra  with  about  the  same  discrepancy  as  for  case   1,      We  also  note 
that  these  computed  values  appear  to  be,    for  most  cases,    shifted  to 
the  right  by  a  small  amount.     No  explanation  other  than  a  general  im- 
precision in  the  determination  of  the  transfer  function  has  been  found 
for  this  shift.     It  is  remarkable  to  note  that  the  rms  values  for  all 
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these  cases  remain  relatively  close  together.     The  following  table  con- 
tains all  rms  values  measured  and  computed  for  all  cases  of  series  I 
and  series  II. 


RMS   VALUES   OF    THE    ACCELERATION   AT    POINT    0   FOR    LOADS 
AT    POINTS    1    AND    2    FOR    SERIES   I   AND    II 


Difference 

r 

Series 

Cases 

Computed 

Measured 

above 

rms  in  g's 

rms  in  g'  s 

measured 

v^ 

% 

1 

3.8 

3.9 

-3 

15.3 

2 

4.45 

4.2 

+  7 

16.3 

I 

3 

4.45 

4.  15 

+  7 

16.  1 

4 

4.  3 

4.  1 

+  5 

15.  1 

5 

4.6 

4.2 

+  10 

16.6 

6 

5.  0 

4.6 

+  9 

18.2 

7 

4.6 

4,2 

+  10 

16.0 

n 

8 

4.55 

4.8 

-6 

19.0 

9 

4.35 

4.  15 

+  5 

16.8 

3.     Effect  of  Correlation 
Special  attention  is  now  given  to  cases  1  and  5.     Case  5  has  a 
completely  uncorrelated  loading  as  opposed  to  case   1,    which  is  com- 
pletely correlated.     In  both  cases  the  pilot  signals     i/Qt)  used  to  drive 
the  two  amplifiers,    nos.    1  and  2,    were  selected  from  pressure  rec- 
ords made  previously  at  the  same  location  of  the  air  jet  [15]   so  that 
they  all  have  the   same  autocorrelation  function.     In  case   1  the  forces 
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fi  (t)  and    b  ^i)  were  obtained  by  feeding  the  same  record      J(t)   simul- 
taneously into  the  power  amplifiers,    nos.    1  and  2.     In  case  5  two  dif- 
ferent records        o  (1}  and         Si^i)  taken  at  two  different  times  were  con- 
sidered so  that  the  cross-correlation  between  them  was  lost,    contrary 
to  case   1  where  the  cross-correlation  was  maximum  since  only  one 
single  record  was  used.      This  idealistic   situation  was  only  partially 
altered  by  the  annplifiers  axid  the  feedback  effect  as  shown  in  Figures 
27  and  35. 

Output  acceleration  spectra  of  cases   1  and  5  are  shown  in 
Figure  37  in  an  attempt  to  point  out  the  effect  of  the  cross-correlation. 
The  distribution  of  frequencies  is  different  for  the  two  cases.      There 
is  more  power  distributed  in  the  higher  frequency  range  for  the  un- 
correlated  load  than  for  the  correlated  case. 

This  difference  has  been  explained  by  analyzing  the  contribu- 
tion of  each  term  in  the  computation  of  the  response  in  case   1.     Corre- 
lated and  uncorrelated  parts  were  isolated  in  the  computation  as  stated 
in  Section  8,    Chapter  I.      These  two  terms  are  represented  in  Figure  38. 
It  can  be  seen  that  above  about  700  cps  the  correlated  part  is  of  the 
same  order  of  magnitude  as  the  uncorrelated  part.      This  shows  clearly 
that  the  cross-correlation  between  the  input  forces  has  a  real  impor- 
tance and  that  neglecting  it  would  lead  to  an  erroneous  spectrum. 
Moreover,    above  700  cps  the  total  spectrum  is  the  difference  of  two 
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Figure  38.     Connputed  Accelerations,    Correlated  and  Uncorrelated 
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large  numbers,    which  explains,    at  least  partially,    the  larger  discrep- 
ancy between  measured  and  connputed  values  in  this  region. 

4.     Series  II  Loading 
In  this  series  a  more  hypothetical  loading  was  investigated. 
The  two  coils  were  hooked  up  in  such  a  way  that  for  a  direct  current 
through  the  coils,    one  coil  would  push  while  the  other  would  pull.     Then 
a  single  random  signal      J  (t)  was  sent  through  amplifier  no.    1  and  the 
same  random  signal,   but  delayed  by  a  tinne     "^  ,    JCt-T:]     ,    was  sent  to 
the  second  amplifier.     Four  cases  of  time  delay,     T   =  0,     T   =   0.  5  mil- 
lisec,     ^   =    1  millisec,     T   =   2  millisec,    were  tried,    leading  to  the 
cross-correlations  indicated  in  Figure  39.      The  rms  values     rj      and    r^ 

•^     of  the  input  forces     p,C^)  a.nd    p„(t)  are  also  indicated  in  Figure  39, 
and  for  sinnplifying  the  computation,    the  autocorrelation  was  assunried 
to  be  the  same  as  for  case   1. 

The  resulting  acceleration  spectra  were  measured  and  com- 
puted with  the  sanne  technique  as  above.      The  results  are   shown  in 
Figures  40  through  43.      The  discrepancies  between  measured  curves 
and  computed  values  are  of  the  same  order  as  for  the  previous  series. 
However,    the  curves  indicate  a  definite  effect  of  the  time  delay   X    on 
the  shape  of  the  spectra.     This  series  shows  that  not  only  the  majci- 
mum  of  the  cross-correlation  is  important,    but  also  that  its  time  posi- 
tion on  the  time   scale  of  the  correlation  function  plays  a  determinant 
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Figure  40.     Power  Spectral  Densities  of  Acceleration,    Case  6 
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Figure  41.     Power  Spectral  Densities  of  Acceleration,    Case  7 
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Figure  42.     Power  Spectral  Densities  of  Acceleration,    Case  8 
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Figure  43.     Power  Spectral  Densities  of  Acceleration,    Case  9 
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role.     It  is  noteworthy  to  mention  here  that  similar  kinds  of  experi- 
ments were  previously  conducted  on  a  string  [3]. 

5.      Conclusion 

In  this  experimental  investigation  it  has  been  shown  for  sev- 
eral loading  configurations  that  once  the  transfer  function  of  a  canti- 
lever beam  has  been  determined,    at  least  for  a  finite  frequency  band, 
then  it  is  possible  to  predict  the  response  to  a  two-point  random  load- 
ing with  reasonable  accuracy. 

Also,    it  has  been  found  that  the  cross-correlation  of  the  load- 
ing plays  a  definite  role  in  the  determination  of  this  response,    and 
should  not  be  ignored  when  the   spectral  distribution  is  of  interest. 
However,    the  root  mean  square  value  of  the  response  has  been  found 
to  be  relatively  insensitive  to  the  cross-correlation  effect. 

These  experiments  have  shown  that  it  is  perfectly  feasible  to 
determine  the  response  of  a  structure  to  a  random  space-time  discrete 
input  by  a  semiexperimental  method,    at  least  for  a  two-point  loading. 
The  extension  to  the  case  of  a  large  number  of  loading  points  does  not 
present  any  mathematical  difficulties,    but  the  amount  of  data  process- 
ing and  numerical  computation  would  become  rather  large.     In  this  re- 
spect an  automatic  correlator  or  even  an  electronic   systenn  for  a  direct 
determination  of  the  cross-power  spectrum  would  be  desirable. 
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